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In the experimental tests of gravity, there have been considerable interests in the pos- 
sibility of intermediate-range gravity. In this paper, we use the earth-satellite mea- 
surement of earth gravity, the lunar orbiter measurement of lunar gravity, and lu- 
nar laser ranging measurement to constrain the intermediate-range gravity from A = 
1.2 X 10^m-3.8x lO^m. The limits for this range are a = 10 * — 5x10 which improve 
previous limits by about one order of magnitude in the range A = 1.2 X 10^ m— 3.8 X 10* m. 

1. Introduction 

To test theories of gravity and to probe its origin, there have been considerable 
interests in the search for non-Newtonian gravityP^Both composition-dependent and 
composition-independent deviations from Newtonian gravity have been searched for. 
Composition-independent deviations from the Newton's inverse square law between 
two point masses M and m can be expressed in terms of a distance-dependent 
gravitational "constant" G{r) as 

G(r) = FgravT^/Mm, (1) 

where r is the separation distance between two masses. Wagonar,^ Fujii, and 
O'HanloiP have provided a theoretical basis for this variation and proposed a — fi 
model as 

G(r) =Ge[l + a(l+M^)e"n, (2) 

where Gc, a, and /i(= A^^) are constants, and G(r) as a function of r is shown in 
Fig. 1. As shown in the figure, G(r) goes to Gc[l + a] with decreasing distance , 
and goes to Gc with increasing distance. 
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Since the proposal of the a — fi model, there have been theoreti cal ana lyses 
and experimental efforts to constrain the variable ranges of a and ^.HEHZI The 
cumulated constraints together with the result obtained in this paper are shown in 
Fig. 2. 

In this paper, we use the satellite measurement of earth gravity, the lunar orbiter 
measurement of lunar gravity, and lunar laser ranging measurement to constrain the 
intermediate-range gravity from A = 1.2 x 10^ m — 3.8 x 10^ m. In section 2, we 
discuss the value of the GM^arth from satellite measurement of earth gravity. In 
section 3, we discuss the value of GMmoon from the lunar orbiter measurement of 
lunar gravity. In section 4, we discuss the G' Mearth+moon obtained from the lunar 
laser ranging experiment. In section 5, we fit the a parameter of the a — fi model 
to these values to constrain the intermediate-range gravity. In section 6, we look 
into possible improvements in the future. 

2. Measurement of GMearth from Satellite Orbit Determination 

The geocentric gravitational coefhcient of the earth [GM earth) is a key parameter of 
earth gravity model. It is determined by observing the influence of this parameter on 
the motion of earth satellites. LAGEOS (LAser GEOdynamics Satellite), launched 
in 1976, was designed to minimize the effects of non-gravitational forces. Because 
of its high altitude (5900 km), the high-Z components of the earth's geopotential 
are greatly attenuated. The result of GM^arth determined by LAGEOS satellite 
has been improved consistentlyEI In 1985, a value of 398600.440 ± 0.002 km^/sec^ 
for GMearth was determined based on 8 years data of laser ranging to LAGEOS 
by Tapley et a This imp rov ed by one order of magnitude than that of GM^arth 
(398600.44 ± 0.02 km^/sec^lEII 

determined using laser ranging to 4 near-earth satel- 
lites in 1978. In 1989, Ries et reported a solution for GM^arth obtained from 
LAGEOS laser ranging and also from a multi-satellite solution to be 398600.4405 ± 
0.001 km^/scc^. 

In 1992, Ries et aL^^ found that the value of the correction for the offset between 
the LAGEOS center-of-mass and the effective reflecting surface should be shifted 
by 1 1 mm. After the center-of-mass offset error had been corrected, they improved 
determination of GM^arth by two methods. The first method is the LAGEOS-only 
solution, where the data set consisted of laser ranges from over 60 stations spanning 
the five-year period from November 1986 to November 1991. The second one is that 
laser range, Doppler and altimeter observations from 17 near-earth satellites are 
combined with surface gravity data in a solution for GMearth- The effects of general 
relativity are taken into account in these two methods. The determined values for 
GMearth by two abovc- mentioned methods are remarkable in agreement with each 
other. The value of GMearth determined there is 398600.4418 ± 0.0008 km^/sec^ in 
SI units or 398600.4415 ± 0.0008 km^/sec^ in TT (Terrestrial Time) units. The scale 
difference of Lq 6.969290134 x 10"^° between the TT and SI units results in a 
change of AGMearth = 0.0003 km^/sec^.Eini The value of GMearth, including the 
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mass of the earth's atmosphere, has been estimated with 1-cr uncertainties in the 
measurement data. The rms of the residuals from the short-arc fit was 2.8 cm.l^ 
UT/CSR (University of Texas/Center of Space Research) has obtained a solution of 
the gravity field complete to degree and order 50 by the second method. Since then 
more endeavors have been taken into developing earth gravity model with higher 
degree and order, while the value of GMgarth has not been improved. In fact, the 
value is still adopted by the IERS(2003) standard released recently. "'^^ 

A new EGM96 gravity Model, which was proposed by NASA (National Aeronau- 
tics and Space Administration) , NIMA (National Imagery and Mapping Agency), 
and OSU (Ohio State University ) in 1993, is developed to support terrestrial and 
extraterrestrial scientific endeavors in connection with its associated global geoid. 
That is one of most accurate earth gravity models and its geoid undulation with 
respect to the WGS84 (World Geodetic System 1984) ellipsoid is about 0.5 m to 
1.0 m. 

CHAMld (CHAUenging Minisatellite Payload), which was launched on July 
15, 2000, was proposed by DLR (Deutschen Zentrum fiir Luft- und Raumfahrt) 
in 1994 and managed by GFZ (GeoForschungsZentrum Potsdam). Determining 
the earth gravity field is one of main scientific objectives of this project. The 
CHAMP project improves by two orders of magnitude in precision for determining 
the earth gravity field. Thus CHAMP obtains low frequency characteristics of high- 
precision earth static gravity field model and their variations with time. EIGEN- 
CHAMPOS^ is a CHAMP-only gravity field model derived from CHAMP GPS 
satellite-to-satellite and accelerometer data out of the period from October 2000 
to June 2003. The accuracy of EIGEN-CHAMP03S is about 5 cm and 0.5 mgal in 
terms of geoid heights and gravity anomalies, respectively, with space resolution at 
A/2 = 400 km. 

GRACE,E1 which was launched on March 17, 2002, is a joint project between 
NASA and DLR. The primary objective of the GRACE mission is to provide with 
unprecedented accuracy estimates of the global high-resolution models of the earth's 
gravity field for a period of up to five years. With the increase of GRACE mea- 
surement data, several Gravity field models derived by GRACE-only data or by 
combined GRACE and CHAMP data have been developed, including GGMOIS 
(complete to degree/order 120), GRACE02S (complete to degree/order 150), and 
GGM02S (complete to degree/order 200). IHl 

The value of GMearth accepted by all of above-mentioned earth gravity models 
[IERS(2003), EGM96, EIGEN-CHAMP03S, GGMOIS, GRACE02S, GGM02S] is 
that one determined by Ries et al^ in 1992, i.e., 

GMearth = 398600.4418 ± 0.0008 km^/sec^. (SI units) (3) 

3. Measurement of GMmoon from the Orbit Determination of Lunar Or- 
biters 
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Like earth, the GMmoon of the Moon is also a key parameter of Moon gravity 
field model. In early period, Lunar Orbiter series and Apollo series space projects 
from later 1960s to early 1970s have obtained large numbers of data for developing 
Moon gravity field model. The Clementine mission was launched on January 25, 
1994. One of its key objectives is to establish a high-resolution gravity field model 
of the Moon by measurements of perturbations in the motion of the spacecraft 
to infer the lunar gravity field. The GLGM-2 (Goddard Lunar Gravity Model 
2 jlSt is the final solution of the Clementine gravity model. The gravity solution is 
based on Doppler tracking of Lunar Orbiters 1 to 5, the Apollo-15 and Apollo-16 
subsatellites, and the Clementine spacecraft .1^^ In order to obtain a more realistic 
estimate of the uncertainty of gravity field, the solution was calibrated by the use of 
subset solutions, as described by Lerch.l^ The value of GMmoon for the GLGM-2 
is 4902.80295 ± 0.00224 km^/sec^. 

On January 6, 1998, NASA launched Lunar Prospector spacecraft, which is the 
third Discovery series mission. For about one- year. Lunar Prospector was placed in 
a near circular orbit at an altitude 100 km and remained in this 2-hour orbit with 
90° inclination. This provided global coverage for the lunar gravity experiment 
every 14 days. On December 19, 1998, the altitude of spacecraft was reduced to 
an average of 40 km to calibrate the gravity field in preparation for an even lower 
extended mission. After January 29, 1999, Lmiar Prospector began its extended 
mission when the spacecraft was lowered to an average of 30 km to obtain higher 
resolution gravity data."^*-* 

Since there is no direct measurement of the lunar farside from Lunar Prospector 
or any other mission, gravity details for the farside of the moon are quite limited. 
However, it has almost no infiuence on the value of GMmoon- The lOOth-degree 
lunar gravity models (LPIOOJ and LPIOOK) extract most of the information from 
the nominal 100 km altitude. In these Lunar Prospector models series, the values 
of GM^oon are hsted in Table 1, where the value of LPIOOK model, 4902.800238 ± 
0.000206 km^/sec^, is the most accurate one, i.e., 

GMmoon = 4902.800238 ± 0.000206 kmVsec^. (SI units) (4) 

Table 1. Key parameters of Lunar Prospector series gravity field modelsl^3 



Model 


G A'Imoon 






Degree / order 




(km-^/sec^) 








LP75D 


4902.801374 




0.00031 


75 


LP75G 


4902.800269 


± 


0.000233 


75 


LPIOOJ 


4902.800476 


± 


0.000209 


100 


LPIOOK 


4902.800238 


± 


0.000206 


100 



4. Determination of G' Mearth-\-moon from the Lunar Laser Ranging Exper- 
iment 
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The accuracy of Lunar Laser Ranging have improved during last 35 years to 2 
cm. The accuracy of Lunar Laser Ranging data from 1970 to 1997 adopted by DE 
405, is about 2-3 cm; the relative accuracy of ranging is better than 8 x 10"^^. The 
latest DE 410 has been created especially for the positions of Mars and Saturn, 
and it includes recent determinations of GM values of the earth and the moon. In 
fact, the GM values discussed in section 2 and 3 are adopted by DE 410. In the 
earth-moon center-of-mass reference frame, the Newtonian equation of the motion 
of moon is 

r = — r, (5) 

where G" is the gravitational constant at earth-moon distance, i.e., 386000 km. 

In DE 405, the GM constant of earth-moon system, which is in TDB (Barycentric 
Dynamical Time) Scale, derived using Lunar Laser Ranging data is 1121 



G' M,arth+-moon = 0.8997011347 X 10-9 AUVday3 

= 403503.23348 W/sec^. (TDB scale) (6) 

When this is converted to SI units, it sh ould b e multiplied by a factor of 1 + Lb with 
Lb = 1.55051976772X 10-8±(2xl0-"),inilland the value is 403503.23974 km^/sec^. 
The uncertainty of this value is not given in the literature of DE 405. However, 
we can estimate the uncertainty conservatively using the accuracy of Lunar Laser 
Ranging data explained at the beginning of this section, i.e., 8 x 10^^^ in fraction of 
r; for ^^q,^P = 2 x the fractional uncertainty of G' Mearth+moon is 1.6 x 10^^^ 
or ±0.00006 km'^/sec^. Putting the value and uncertainty together, G' Mearth+jnoon 
is 

G' Mearth+moon = 403503.23974 ± 0.00006 kmVsec^. (SI units) (7) 
5. Constraint on the a — jj. model 

Comparing with the separation distance from earth to moon, the distance of (1.8 — 
12) X 10^ m, on which the gravity field models based on earth's satellite or moon's 
satellite are constructed, is rather short. For these short distances, we can obtain 
approximatively 

Mraoon/M^arth ^ 0.0123000371 ± (5 X 10-1"), (8) 

from (3), (4) in section 2 and 3. DE 410 also gives this mass ratio as 0.0123000371; 
the same as (8) without quoting an uncertainty. However, DE 405 adopts the value 
as 0.0123000383. 

From (7) and (8), we obtain 
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G'Mearth = 398600.4395 ± 0.00022 W/sec^. (9) 

In (9), G' is measured at the separation distance from earth to moon of r„j = 
3.86 X 10^ m. According to section 2, at the distance of earth satelhte of — 
1.2 X 10^ m, we have 

GMearth = 398600.4418 ± 0.0008 W/sec^ (10) 
Dividing (10) by (9), we obtain 

G/G" 1.000000006(±6 X 10"^). (3 - cr) (11) 

According to a — fi model, now we use (11) to constrain the parameter space of 
intermediate range gravity (2). Assuming < A < and setting 



(11) becomes 



M/(Ai,r)-(l + He-'^^ (a^ = A-1) (12) 

^-f^f^-lM6±6)xlO^ (13) 
In (13), a is small and (13) can be approximated as 

[■^ifi, r,) - ^'(/., rm)]a = (6 ± 6) x 10'^. (14) 

Thus, 

^ 12 X lO-'^ 

[vI/(M,r,)-*(^,r™)]- ^ > 

According to (15), for a given A in the range 1.2 x 10^ m — 3.8 x 10* m, the a 
parameter of the a — fi model is constrained as in Fig. 3. This new result improves 
the constraint in the range 1.2 x 10^ m — 3.8 x 10* m on Fig. 2 by about one order 
of magnitude. 

6. Discussions and Outlook 

In this paper, we have used the satellite measurement of earth gravity, the lunar 
orbiter measurement of lunar gravity, and laser ranging measurement to constrain 
the intermediate-range gravity from A = 1.2 x 10^ m — 3.8 x 10* m. For the satel- 
lite measurement of earth gravity, the geodesy mission GOCE (Gravity field and 
steady-state Ocean Circulation Explorer^ is planned to be launched in 2006. The 
scientific objective is the mapping of the Earth's static gravity field with a very high 
resolution and accuracy on a global scale. GOCE is a drag-free mission, flown in 
a circular and sun-synchronous orbit at an altitude between 240 km and 250 km. 
A successful mission presupposed, GOCE would finally deliver the Earth's gravity 
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field with a resolution of about 70 km half-wavelength and a global geoid with an 
accuracy of about 1 cm. When this and other geodesy missions are implemented, 
the GMearth wiU be determined more accurately. 

SMART-lpSI the Europe's first lunar mission launched in 2003, is playing a 
vital role in developing cutting edge technologies that will be a major part of the 
future of lunar and planetary science. Japan plans to launch LUNAR-A mission in 
2005,12^ whose scientific objective is to explore the lunar interior using seismometry 
and heat-flow measurements. However, improvement in lunar gravity model is not 
a goal of SMART-1 and LUNAR-A. 

Japan is planning to launch SELENE (SELenological and ENgineering Explorer) 
using H-2A rocket in 2006.123 The project will use a main orbiting satellite and two 
subsatellites to obtain highly accurate lunar gravity model. After that, it will be 
possible to construct a global gravity model of the Moon using spherical harmonics, 
and the value of GMmoon will be improved. 

Lunar laser ranging at CERGA is under renovation and improvement.!^ APOLLO 
program will use the 3.5 m telescope at the Apache Point Observatory for lunar 
laser ranging.l2HI These experiments will push Lunar Laser Ranging into millimeter- 
accuracy range. As a result, constrains in this paper will also be improved. 

For small A, analysis of satellite ranging data together with GOCE will improve 
the constraints. For lager A, we look into ASTROD I and ASTROD missions 
which are currently under pre-phase A study. ASTROD I and ASTROD are to 
implement deep-space pulse and interferometric laser ranging in the solar system, 
to test relativistic gravity and to increase the sensitivity of solar, planetary and 
asteroid parameter determination by 1 to 5 order of m agnitudes . This will improve 
constraints on a for larger A significantlyP^^^^ESSni 

In this paper, we have used data from earth gravity, lunar gravity, and lunar laser 
ranging separately to constrain the strength of intermediate-range gravity. In the 
future, we will use the a — ^ model as the underlying gravity theory for a grand fit of 
all three kinds of data together. With increased accuracy of gravity measurements 
and lunar laser ranging, we expect the constraints on intermediate-ranging gravity 
to be significantly improved. 
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